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The quantum Jarzynski equality is an important theorem of modern quantum thermodynamics. 
We show that the Jarzynski equality readily generalizes to relativistic quantum mechanics described 
by the Dirac equation. After establishing the conceptual framework we solve a pedagogical, yet 
experimentally relevant, system analytically. As a main result we obtain the exact quantum work 
distributions for charged particles traveling through a time-dependent vector potential evolving 
under Schrodinger as well as under Dirac dynamics, and for which the Jarzynski equality is verified. 
Special emphasis is put on the conceptual and technical subtleties arising from relativistic quantum 
mechanics. 
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I. INTRODUCTION 

The Jarzynski equality [1] together with subsequent 
Nonequilibrium Work Theorems, such as the Crooks fluc¬ 
tuation theorem [2], are undoubtedly among the most 
important breakthroughs in modern Statistical Physics 
[3]. In traditional thermodynamics the only processes 
that are fully describable are infinitely slow - equilibrium 
- processes [4]. For all realistic, finite-time - nonequi¬ 
librium - processes the second law of thermodynamics 
constitutes an inequality, only stating that some portion 
of the entropy is irreversibly dissipated into the environ¬ 
ment. Jarzynski showed that for isothermal processes the 
second law of thermodynamics can be formulated as an 
equality, no matter how far from equilibrium the system 
is driven [1], (exp (—/3kF)) = exp {—/3AF). Here (3 is the 
inverse temperature of the environment, and AF is the 
free energy difference, i.e., the work performed during an 
infinitely slow process. The angular brackets denote an 
average over an ensemble of finite-time realizations of the 
process characterized by their nonequilibrium work W. 

The discovery of these so-called fluctuation theorems 
effectively opened a new field of contemporary research 
[5, 6]. For small, but classical systems the Jarzynski 
equality is a universally valid theorem [7] , which has been 
experimentally verified in a variety of systems [8-11]. For 
quantum systems, however, the situation is more compli¬ 
cated. The major conceptual obstacle is how to gener¬ 
alize the classical notion of thermodynamic work to the 
quantum domain. In particular, quantum work is not an 
observable in the usual sense, as there is no hermitian op¬ 
erator, whose eigenvalues are given by the classical work 
values [12-16]. 

For isolated quantum systems evolving under unitary 
dynamics the so-called two-time energy measurement ap¬ 
proach has proven to be practical and powerful. In this 
paradigm, quantum work is determined by projective en¬ 
ergy measurements at the beginning and the end of a 
process induced by an externally controlled Hamiltonian. 
Although this approach has been verified experimentally 
[17-21] and has led to the development of thermodynamic 


quantum devices [22-24] , the paradigm cannot be consid¬ 
ered entirely satisfactory as it relies on a rather invasive 
procedure - projective measurements - and is restricted 
to isolated systems. 

Thus, modern quantum thermodynamics has been at¬ 
tempting to overcome these restrictions: On the one 
hand, researchers have generalized the two-time energy 
measurement approach to less invasive procedures such 
as generalized measurements [25-30] , or to quantum sys¬ 
tems that are less “isolated” such as in T^T-symmetric 
quantum mechanics [31]. On the other hand, various no¬ 
tions of quantum work and entropy production for gen¬ 
eral, open quantum systems have been proposed [32-37], 
which however all lack the desired universality of notions 
from traditional thermodynamics. 

Nevertheless, due to its simplicity and practicality for 
isolated quantum systems a great deal of research has 
been dedicated to a careful study of the quantum work 
statistics from two-time energy measurements. For in¬ 
stance, the quantum work distribution has been com¬ 
puted for time-dependent oscillators [38-40], a particle 
in a time-dependent box [41] , quantum Ising chains [42- 
45], the Landau-Zener model [46], noninteracting bosons 
and fermions [47], diatomic molecules [48], etc. 

However, to the best of our knowledge all prior work 
has focused on non-relativistic quantum systems, while a 
generalization of the Jarzynski equality to relativistic en¬ 
ergies has only been proposed for classical systems [49]. 
The present paper aims at closing this gap and reports 
the generalization of the quantum Jarzynski equality to 
particles evolving under the time-dependent Dirac equa¬ 
tion. We will see that the validity of the Jarzynski equal¬ 
ity together with the two-time energy measurement ap¬ 
proach follows directly from the unitarity of Dirac dy¬ 
namics - the only essential requirement [25]. Therefore, 
after briefly establishing the conceptual building blocks, 
we will focus on a pedagogical and illustrative case study, 
namely charged spin-1/2 particles traveling through a 
time-dependent vector potential. 

The purpose of the present study is two-fold: We 
will show that the quantum Jarzynski equality naturally 


2 


holds for dynamics described by the Dirac equation. The 
main part of the following discussion, however, will pro¬ 
vide a “recipe” of how to compute the relativistic quan¬ 
tum work density. Our analysis will put emphasis on the 
technical and conceptual subtleties arising from Dirac’s 
equation, and we will compare our relativistic results 
with the analogous Schrddinger dynamics. 


II. RELATIVISTIC QUANTUM WORK 

We begin by briefly reviewing the paradigm of the 
two-time energy measurement approach, and establish 
notions and notations. Consider an isolated quantum 
system with time-dependent Schrodinger equation 

i/i = i?t IV't) , (1) 

where the dot denotes a derivative with respect to time. 
We are interested in describing thermodynamic processes 
that are induced by varying an external control param¬ 
eter Xt during time r, with Ht = H{Xt). Within the 
paradigm of two-time energy measurements quantum 
work is determined by the following, experimentally mo¬ 
tivated protocol: After preparation of the initial state 
Po a projective energy measurement is performed; then 
the system is allowed to evolve under the time-dependent 
Schrodinger equation (1), before a second projective en¬ 
ergy measurement is performed at t = t. Thus, for a 
single realization of this protocol the quantum work is 
given by 


where |no) is the initial eigenstate with eigenenergy 
e(no, Aq) and |nr) with e{nr, Xr) describes the final eigen¬ 
state. 

The quantum work density is then given by an 
average over an ensemble of realizations, V(W) = 
{S {W — Wno^n^)), which can be rewritten as [25, 34] 

V{W) =^SiW- p K ^ n ^). (3) 


Po = eyip {—l3Ho)/Zo, we have the quantum Jarzynski 
equality [12-15], 

(exp (-/3W)) = exp (-/3AF), (5) 

where AF = F^. — Fq and Ft = —(1//3) In (Zt). 

It is worth emphasizing that the validity of the quan¬ 
tum Jarzynski equality is not restricted to Schrodinger 
dynamics. Rather, it has been shown that Eq. (5) holds 
for all quantum systems, whose dynamics is at least uni- 
tal [25, 31, 50-52]. Unital dynamics preserves the identity 
and can be written as a superposition of unitary quantum 
maps [53]. 

Therefore, to check whether the quantum Jarzynski 
equality holds for Dirac dynamics, one only has to ver¬ 
ify that the corresponding evolution equation describes 
unital dynamics. 

Relativistic quantum mechanics: Dirac equation The 
Dirac equation is a relativistic wave equation, which de¬ 
scribes massive spin- 1/2 particles, such as electrons and 
quarks. In its original formulation for free particles the 
Dirac equation reads [54] 

ih^{p, t) = {coL ■ p + aomc^) ^{p,t). ( 6 ) 

Here, 4' (p, t) is the wave function of an electron with 
rest mass m and momentum p — {pi,P 2 ,Pz), and c is 
the speed of light. In covariant form the matrices a. = 
(Q:i,a 2 ,Qf 3 ) and Oq can be expressed as [55], 

= 7 ° and 7 ° = 7 * . (7) 


The 7 -matrices are commonly expressed in terms of 2 x 2 
sub-matrices with the Pauli-matrices CFx,cry,az and the 
identity I 2 as. 



( 8 ) 


In the latter equation the symbol ^ denotes a sum over 
the discrete part of the eigenvalues spectrum and an in¬ 
tegral over the continuous part. 

To compute V{W) (3) explicitly, one has to determine 
the transition probabilities p (uq —>■ Ur) first. These can 
be written as [25, 31], 

p (no Ur) = tr {n„^ Ur n^o Pon„„u/}, (4) 

where po is the initial density operator of the system 
and Ur is the unitary time evolution operator, Ur = 
T>exp(—i/h dt Fit). Finally, n„ denotes the projec¬ 
tor into the space spanned by the nth eigenstate, which 
becomes for non-degenerate spectra n„ = |n) (n|. 

It is then a simple exercise to show that from the 
definition of 'P{W) (3) and for an initial Gibbs state. 


It is then easy to see that the right side of Eq. ( 6 ), i.e., the 
operator ca -p + aomc^, is hermitian, and consequently 
4 '(p, t) evolves under unitary dynamics. 

Therefore, the quantum Jarzynski equality (5) also 
holds for particles evolving under Dirac dynamics ( 6 ). 
However, we expect the work density function (3) to be 
dramatically different: In contrast to the Schrodinger 
equation (1) the Dirac wave function 4 '(p, t) is a bispinor, 
which can be interpreted as a superposition of a spin-up 
electron, a spin-down electron, a spin-up positron, and 
a spin-down positron [55]. In addition, the momentum 
of Dirac particles is confined by the light cone, whereas 
Schrodinger particles can travel with arbitrary velocities. 

In the remainder of this study we will analyze the con¬ 
sequences of relativistic effects on the quantum work den¬ 
sity for a simple, yet elucidating example. 
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III. CHARGED PARTICLES IN A first order. Thus, a solution of Eq. (10) can be written 

TIME-DEPENDENT VECTOR FIELD as 


For the sake of simplicity we now restrict ourselves 
to a one-dimensional system in cc-direction. In this case 
the 4-component Dirac spinor can be separated into two 
identical 2-component bispinors, which evolve under [56] , 

ih 4'(p, t) = {cp + mc^ a^) 4'(p, t), (9) 

with Pa; = p. We further assume that the system is driven 
by a time-dependent, but spatially homogeneous vector 
potential At. For oscillating At this situation has been 
recently solved analytically [56]. Moreover, Eq. (9) de¬ 
scribes particle-antiparticle production in counterpropa- 
gating laser light, which has been proposed to be observ¬ 
able in an experiment [56]. 

Note that the Dirac equation (9) as any electromagen- 
tic theory is gauge invariant. Here, “gauge invariance” 
means that a whole class of scalar and vector potentials, 
related by so-called gauge transformations, describes the 
same physical situation. In particular, the dynamics of 
the electromagnetic fields and the dynamics of a charged 
system in an electromagnetic background do not depend 
on the choice of the gauge. In the present context this 
means that the energy eigenvalues (II) and (20) do de¬ 
pend on the gauge, whereas the Jarzysnki equality is 
gauge invariant [15]. 


■)p{p,t) = exp y dt' ^ V'feO), (14) 

which follows from inspection. 

Notice that in the case of Schrodinger dynamics the ef¬ 
fect of At manifests itself exclusively as a time-dependent 
phase (14). Shortly, we will see that for the correspond¬ 
ing Dirac equation the situation is much more involved. 

The instantaneous eigenenergies (11) together with the 
eigenstates (12) and the time-dependent solution (14) are 
all ingredients necessary to compute the quantum work 
density (3). In particular, the transition probabilities (4) 
become. 


PS (tTo -t TTt) 


J dp(j){p, TTr)'<p{pA) 


2 

PoK), 


(15) 


with the initial state 

Po(’^o) = exp(-/3es(7ro, Ao))/^^, (16) 

and partition function Zq = f dyro exp (—,0 es( 7 ro, H.o)) 
and ipiPyO) = (/)(p, ttq). Substituting Eqs. (11) and (15) 
into Eq. (3) we finally obtain after a few lines of simple 
algebra. 


A. Schrodinger dynamics 

To build intuition and as a point of reference we treat 
the non-relativistic problem, first. In this case the dy¬ 
namics is described by the time-dependent Schrodinger 
equation, which reads in momentum representation 

ihip(p,t) = (p +At/cf 7p(p,t), (10) 

where At is the vector potential, and we work in units 
for which the elementary charge is set to one, e = 1 . 
Accordingly, the instantaneous eigenenergies are, 

^s(At,At) = ^ {TTt +At/cf ( 11 ) 

2m 

with the corresponding eigenstates, 

(l){p,Trt)=S{p-{Trt + At/c)). (12) 

Here and in the following nt denotes the quantum num¬ 
ber, which is in the present case reduces to the eigen- 
momentum. Note that the eigenstates (12) form an or¬ 
thonormal basis, since 

J dp(j){p,Tri)(j){p,Tr2) = 5 {tti - Tr2) . (13) 

In the present case the time-dependent Schrodinger equa¬ 
tion ( 10 ) reduces to an ordinary differential equation of 


VsiW) = 

^ ( P[2mc^W-{Ao-2ArfY \ 

a/^ |Ao - 2A^| ^ ^ y 8mc2 (Ao - 2At-)2 j 

(17) 

Equation (17) constitutes our first main result. The 
quantum work distribution for charged Schrodinger par¬ 
ticles traveling trough a time-dependent vector potential. 
At, is a Gaussian, which is fully determined by the ini¬ 
tial and final value of At. In particular, Vs{W) is inde¬ 
pendent of the specihc protocol, as At merely induces a 
time-dependent phase (14). As a point of reference and 
for comparison with the Dirac case in the following sub¬ 
section, we plot Eq. (17) in Fig. 1 for low, intermediate, 
and high temperatures. 


B. Dirac dynamics 

In complete analogy to the preceding Schrodinger case 
we now compute the quantum work density for charged 
particles evolving under the time-dependent Dirac equa¬ 
tion, 

ih'^{pA) =[{cp +At) a^] 4'(p,t). (18) 

Equation (18) can be separated into two evolution equa¬ 
tions for the components of the bispinor, 'I'(p, t) = 
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FIG. 1. (color online) Quantum work density, VsiW), for 
charged Schrodinger particles (17) with = 0, Ar = 1, 
m = 1, c = 1, and /I = 10 (blue, solid line), /3 = 1 (purple, 
dashed line), and /3 = 0.1 (red, dotted line). 


('I'i(p, i), 4'2 (p, i)), and we have 

= - {cp + At)"^ + {mc^)^ + ihAt 4'i(p,<), 
4'2(p,t)= ih'^iipA) - {cp +At)'^i{pA) lm(? . 

(19) 


In contrast to the previous case (10) the solution of the 
time-dependent Dirac equation (18) is determined by 
an ordinary differential equation of second order (19). 
Therefore, to find analytical solutions of Eq. (18) we have 
to resort to particular parameterizations of At- For os¬ 
cillating protocols Eq. (18) has been solved in Ref. [56], 
and we will see two further examples in the following. 

Before we turn to specific parameterizations, how¬ 
ever, we note the instantaneous (positive) eigenenergies 
ofEq. (18), 

eoi^t^At) = \J {cTTt + Atf + {mc^f ( 20 ) 


fluctuation theorem can be derived. However, such a 
general theorem only reduces to a generalized Jarzynski 
equality for thermodynamically well-defined situations 
[34]. Nevertheless, for the sake of completeness, antipar¬ 
ticle energies and states can be found in Appendix A. 

We also would like to emphasize that our analysis does 
not neglect the existence of antiparticles completely. We 
merely assume that the initial state is a thermal wave 
packet of particles. The dynamics, however, is described 
by the time-dependent Dirac equation (18), and hence 
governed by both, positive and negative eigenenergies. 
Hence, in particular the transition probabilities (4) are 
governed by both, particle and antiparticle solution. 

a. Linear protocol As a first example we consider a 
linearly parameterized vector potential. 


At = at, 

for which a solution of Eq. (19) is given by 

{i -I- l)(cp + at) 


4'i(p,t) = C'i(p) 

+ C2{p) Dy_i 


y/ah 
{i - l){cp + at) 


( 22 ) 


(23) 


Here, D^{-) denotes the parabolic cylinder function [58] 
of order = iim?c^/2ah, and Ci{p) and C 2 {p) are time- 
independent functions of momentum determined by the 
initial state. 

As mentioned earlier, for Dirac dynamics the solution 
is mathematically more involved, and also the amplitude 
of the wave function depends on the specific parameteri¬ 
zation of At - not only the phase as in the previous, non- 
relativistic case (14). This can be understood as dynam¬ 
ical interference of the two components of the bispinor. 
Nevertheless, the transition probabilities can still be writ¬ 
ten as 


pD (tto 


dp$(p,7r,-) • 4'(p,t) 


2 

p^(7ro), (24) 


where the initial distribution now reads 


and the corresponding, orthonormal eigenstates. 


$(p,7rt) 


{iTt + At/c)) / 1 

yi + (Vif+i^ny ly/nm 



( 21 ) 

where we introduced the notation lit = (ctt^ -|- At)lm(?. 
One easily convinces oneself that these eigenstates, 
^{p,Trt), fulfill the orthonormality condition (13). 

For the following analysis we will assume that in the 
initial state, po, merely particles are present, but no an- 
tiparticles. This assumption is in full agreement with 
typical situations in nature and the mathematical treat¬ 
ment simplifies significantly [57]. We emphasize that this 
assumption merely circumvents the conceptual issue of 
having to define a free energy for antiparticles. It has 
been shown that for any normalized initial state [34] a 


Po (t^o) = exp (-^ e£,(7ro, Ao))/Z^ , (25) 

with which we can compute the quantum work distribu¬ 
tion Pd(W) (3). 

In Fig. 2 we plot the thermal momentum distribution 
for Schrodinger particles (16) together with the distri¬ 
bution for Dirac particles (25). In the Schrodinger case 
we have the well-known (Gaussian) Maxwell-Boltzmann 
distribution. The momentum distribution for relativis¬ 
tic Dirac particles is broader due to the relativistic en¬ 
ergy, mc^, and was first described for classical mechanics 
by Jiittner [59]. The so-called Maxwell-Jiittner distribu¬ 
tion converges towards the Maxwell-Boltzmann distribu¬ 
tion (16) for low temperatures, and decays slower than 
a Gaussian at high temperatures [59]. The major limi¬ 
tation is that the Maxwell-Jiittner distribution neglects 
antiparticles, which however serves our present purpose. 
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Po(?ro) 



^o(W) 



FIG. 2. (color online) Thermal momentum distribution for 
Schrodinger (16) (green, solid line) and Dirac particles (25) 
(dark orange, dashed line) with m = 1, c = 1, Aq = 0, and 
P = 0 . 1 . 



FIG. 3. (color online) Quantum work density, VniW), for 
charged Dirac particles and the linear protocol (22) with a = 
1, T = 1, = 1, m = 1, c = 1, and /3 = 10 (blue, solid line), 

P = 1 (purple, dashed line), and /3 = 0.1 (red, dotted line). 


The resulting work distribution, Vd{W), is plotted in 
Fig. 3, for the same parameters and color coding as for 
Schrodinger particles in Fig. 1. As anticipated the rela¬ 
tivistic effects change the characteristics of the quantum 
work distribution significantly. The most striking differ¬ 
ence with the Schrodinger case in Fig. 1 is that Vd{W) 
has a finite support. This, however, can be understood 
intuitively: Large fluctuations in W are accompanied by 
a large change in momentum. However, the momentum 
is limited by the light cone, and, hence, large fluctuations 
in W are also “cut-off” by the light cone. 

Quantitatively, the finite support can be determined 
by inspecting the transition probabilities (24). It is 
easy to see that pD (ttq — tt^-) cx 6 (ttq — (tt^- -|- At/c)). 
Hence, the only work values contributing to Vd{W) are 
W = eoi'n'o + AtIc.A^') — e£i(7ro,Ao), and therefore 
W € (—2 ar, 2 ar). 


FIG. 4. (color online) Quantum work density, 'Pd(IF), for 
charged Dirac particles and the exponential protocol (26) with 
a = e/(e — 1), r = 1, = 1, m = 1, c = 1, and /3 = 10 (blue, 

solid line), /3 = 1 (purple, dashed line), and /3 = 0.1 (red, 
dotted line). 

Qualitatively, one can understand Fig. 3 by starting 
with the distribution for Schrodinger particles in Fig. 1, 
and “compressing” the curves into the allowed support. 
For low temperatures (blue curve) the left flank is rather 
unaffected, as the distribution lives “far away” from the 
light cone, whereas the right flank is only slightly de¬ 
formed. For higher temperatures the effect becomes more 
prominent, and the distribution becomes “jammed” at 
the edges of the support, i.e., at the light cone. 

b. Exponential protocol To conclude the analysis we 
also compute the quantum work density (3) for a non¬ 
linear parameterization, 

At = a(l - exp(-</r)) . (26) 

Also in this case the time-dependent Dirac equation (19) 
can be solved analytically. However, the solution can no 
longer be written in compact form, and can be found 
in Appendix B. The transition probabilities (24) and the 
initial distribution remain the same by replacing Eq. (23) 
with the expression (B2) everywhere. 

Figure 4 illustrates the resulting quantum work dis¬ 
tributions. We observe that the work distributions re¬ 
sulting from the linear protocol (22) and the exponential 
protocol (26) are nearly indistinguishable - despite the 
solutions (23) and (B2) being complicated expressions of 
special functions. Thus, we conclude that the effect of 
the light cone on the work distribution is more promi¬ 
nent than the interference of the two components of the 
bispinor [60]. 

C. Jarzynski equality 

The validity of the quantum Jarzynski equality follows 
from the unitarity of Dirac dynamics. Nevertheless, it is 
worthwhile to numerically verify its predictions. To this 
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13=10 

P = 1 

P = 0.1 

linear 

exponential 

0.99 

0.96 

0.99 

0.98 

0.99 

1.00 


TABLE I. Numerical verification of the quantum Jarzynski 
equality (27) for the quantum work distributions for Dirac 
particles in Figs. 3 and 4. 

end, we numerically integrated the average exponentiated 
work for the distributions in Figs. 3 and 4. Here, the 
Jarzynski equality becomes, 

(exp {-j5W)) = J dWTaiW) exp (-/3bF) = 1, (27) 

as the free energy difference vanishes. In Table I we sum¬ 
marize the numerical results. We see that the quantum 
Jarzynski equality (27) is, indeed, verified to very high 
accuracy. 

The validity of the quantum Jarzynski equality (27) 
explains another important feature of 7 ^_d(IT). For the 
linear (22) as well as for the exponential (26) protocol 
left and right flank of the distribution are “exponentially 
asymmetric”. This asymmetry constitutes a necessary 
charactertistic of Vd{W) for Eq. (27) to hold. Note also 
that the asymmetry of VoiW) is not an artifact of assum¬ 
ing that the initial state is comprised of only particles, 
but no antiparticles. We emphasize again that the exis¬ 
tence of antiparticles is implicit in our treatment as the 
dynamics is described by Eq. (18). 


[56]. Nevertheless, for all experiments for which Eq. (9) is 
valid our results could be readily verified, where one only 
would have to additionally measure the momentum dis¬ 
tribution. From the momentum distribution one would 
compute the transition probabilities (24), and build the 
quantum work distribution (3) from a histogram. This 
procedure is fully analogous to the cold ion trap ex¬ 
periment, that verified the quantum Jarzynski equality 
[17, 21]. 


B. Summary and outlook 

Our present analysis has extended the scope of quan¬ 
tum stochastic thermodynamics to relativistic energies. 
We have shown that not only does the quantum Jarzynski 
equality hold for Dirac dynamics, but we also have pro¬ 
vided a step-by-step “recipe” of how to compute the rela¬ 
tivistic work distribution. For the sake of clarity and due 
to its mathematical simplicity we focused on free, charged 
particles traveling through a time-dependent vector po¬ 
tential. Another recent reference proposed to study pair 
production in a slightly more complicated, but also more 
realistic system including a scalar potential [61]. Our 
analysis could be straightforwardly applied to the situa¬ 
tion of Ref. [61] under the expense of having to compute 
VniW) fully numerically. 
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Appendix A: Antiparticle energy and eigenstate 

The instantaneous antiparticle solution of the time- 
dependent Dirac equation (18) is given by 


A. Experimental relevance 

Only recently, Fillion-Gourdeau et al. [56] studied the 
same model system in the context of pair production in 
counterpropagating laser light. However, Ref. [56] not 
only solves Eq. (18) analytically for an oscillating pa¬ 
rameterization of At, but also provides relevant values 
for the field strength, for which the dynamics could be 
observed in an experiment. It is worth emphasizing that 
Eq. (9) is only an approximate description of the real 
physical situation with a clearly defined range of validity 


$“(T,7rt) 


^{p- {-Kt + At/c)) 

+ ("s/n? + 1 + Ht^ 


with eigenenergies 


-1 


•^/Hj + 1 + H/ 


(Al) 


e^(7rt. At) = (cTTt J- Atf + {mc^f . (A2) 

Note that the eigenenergies for antiparticles are negative. 
Thus, the Maxwell-Jiittner distribution (25) is ill-defined 
for positive temperatures. 
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Appendix B: Analytical solution of time-dependent Dirac equation 

For the exponential protocol, 

At =a(l-exp(-t/r)) 

a solution of the time-dependent Dirac equation (19) is given by 


(Bl) 


{p, t) = exp ^ (t a/ {mc^Y + {cpY + ‘2-pca + + ar 

X C*i(p) U (^cp + a — r-y/ {mc^)^ + {cpY -\- 2cpa + a'^'j , 1 -I- -I- (cp)^ -I- 2cpa + 

+ C 2 ip) L (^cp -I- a — r-y/ {mc^Y + {cpY + 2cpa + + {cpY + 2cpa + 




(B2) 


Here, U{-, •, •) is the Kummer function, and L(-, •, •) denotes the Laguerre polynomial 
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